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ABSTRACT 

In order to assess qualitatively the ejecta geometry of stripped-envelope core-collapse 
supernovae, we investigate 98 late-time spectra of 39 objects, many of them previously 
unpublished. We perform a Gauss-fitting of the [O i] AA6300, 6364 feature in all spectra, 
with the position, full width at half maximum (FWHM) and intensity of the A6300 
Gaussian as free parameters, and the A6364 Gaussian added appropriately to account 
for the doublet nature of the [O i] feature. On the basis of the best-fit parameters, the 
objects are organised into morphological classes, and we conclude that at least half of 
all Type Ib/c supernovae must be aspherical. Bipolar jet-models do not seem to be 
universally applicable, as we find too few symmetric double-peaked [Oi] profiles. In 
some objects the [O i] line exhibits a variety of shifted secondary peaks or shoulders, 
interpreted as blobs of matter ejected at high velocity and possibly accompanied by 
neutron-star kicks to assure momentum conservation. At phases earlier than ~ 200 d, a 
systematic blueshift of the [O i] AA6300, 6364 line centroids can be discerned. Residual 
opacity provides the most convincing explanation of this phenomenon, photons emitted 
on the rear side of the SN being scattered or absorbed on their way through the ejecta. 
Once modified to account for the doublet nature of the oxygen feature, the profile of 
Mgi] A4571 at sufficiently late phases generally resembles that of [Oi] AA6300,6364, 
suggesting negligible contamination from other lines and confirming that O and Mg 
arc similarly distributed within the ejecta. 
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1 INTRODUCTION 

The geometry of stripped-envelope core-collapse supernova 
(CC-SN) ejecta has been scrutinised for about ten years, 
since the association of S N 1998bw with the nearby j- 
ray burst GRB 980425 (|Galama et al.l ll99Sl) . Together 
with subsequent ex amples of SN-GRB associations (see e.g. 
iFerrero et aH l2006 for an overview) this suggested that at 
least some Type Ib/c supernovae (SNe Ib/c) may be pow- 
ered by the same engine as long-duration GRBs, and thus 
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that their ejecta may show large-scale asphericity along 
an axis defined by the GRB-jet. Indeed, the nebular spec- 
tra of SN 1998bw exhibited properti es which could not 
be explained with spherical symmetry (jMazzali et al.ll200ll ; 
iMaeda et aI1l2002l l. Instead, a model with high- velocity Fe- 
rich material ejected along the jet axis, and lower-velocity 
O in a torus perpendicular to this axis, was proposed. From 
this geometry a strong vie wing-angle depend ence of nebular 
line profiles was obtained ()Maeda et al.ll2002l ). Of particular 
note to this study, this model suggests that double-peaked 
O lines should be observed if viewed from a direction per- 
pendicular to the jet. 
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Such a profile was first observed in SN 2003jd 
l|Mazzali et alj 120051 ). an important step towards a coher- 
ent picture. However, whether large-scale asphericity is 
found only in SNe 1c connecte d with a GRB (probably 
only a few percent of all SN e Ic; IPodsiadlowski et al.l 12004 ; 
iGuetta fc Delia Valid 120071 ) . and how large the degree of 
asphericity in 'norm al' SNe Ib/c actually is, are debated. 
iMaeda et al.l l|2008l ) recently studied a sample of nebular 
spectra of 18 stripped-envelope CC-SNe and found a large 
fraction of double-peaked [O i] AA6300, 6364 line profiles, 
consistent with about half of all SNe Ib/c bei ng strongly, or 
all of them moderately, aspherical. Similarly. iModiaz et al.l 
(2008) analysed late-time spectra of 8 stripped CC-SNe, con- 
cluding that asphericity is ubiquitous in all these events, 
not only the hyperenergetic ones. It should be noted that 
the Maeda et al. and Modjaz et al. works focus on devia- 
tions from sphericity on global scales, as opposed to small- 
scale dumpiness of the ejecta that results in fine-structured 
emission lines and requires fa irly high- resolution nebular 
spectra to be studied (see e.g. Filippenko fc Sargentl 1 1983 ; 
ISpvromilioll 19941 ; iMatheson et al.ll200d ). 

In this work we c onduct a study sim ilar to that of 
IMaeda et ail (|2008l ) and lModiaz et all |2008l ), but based on 
a larger sample of SNe, considering virtually all nebular 
SN Ib/c spectra we could access from the literature, com- 
plemented by 26 previously unpublished spectra from the 
Asiago Supernova Archive and recent observations carried 
out at the ESO Very Large Telescope (VLT). The work is 
organised as follows: in Section [5] we present the entire SN 
sample, discuss the selection criteria for spectra to be in- 
cluded, and discuss in more detail those spectra which were 
previously unpublished. Section [3] concentrates on the link 
between the ejecta geometry and observed line profiles, mo- 
tivates the choice to focus on [O i] AA6300, 6364, addresses 
the complications arising from its doublet nature, and in- 
troduces the fitting procedure employed to gain qualitative 
insight into the ejecta morphology. The results of this fitting 
are analysed in Sections and trying to find an explana- 
tion for the mean blueshifts of the line's centroids at phases 
< 200 d, and dividing the objects into different classes on 
the basis of the best-fit parameters. Individual objects with 
interesting line profiles are discussed more deeply in Sec- 
tion HJ1 while Section [7] extends the analysis to the profile of 
Mgl] A4571 and its comparison to that of [Oi]. Finally, a 
brief summary of the main results is given in Section [8] 



2 THE SAMPLE OF SN Ib/c SPECTRA 

Our goal is to compare a large set of late-time spectra of 
stripped-envelope CC-SNe, concentrating on what can be 
learned about the ejecta geometry by studying the profiles 
of nebular emission lines, in particular [O i] AA6300, 6364 
(the motivation to focus on this line is given in Section [3]). 

Given the statistical approach of this study, we use a 
simple fitting procedure (for details see Section [3} . Com- 
pared to full spectral modelling this method has the ad- 
vantage of being fast, capable of dealing with complex pro- 
files, and independent of an accurate flux calibration of the 
spectra, thus allowing it to be applied to a large number of 
spectra. 

Nebular emission features in SNe Ib/c typically start 



to emerge about two months past maximum light, but at 
that epoch the SN flux is still dominated by photospheric 
emission. For this reason we included in our sample only 
spectra which were taken more than ~ 90 d after maximum 
light, which, assuming typical rise times, corresponds to 100 
or 110 d after explosion. At those phases there still is an 
underlying photospheric continuum, but this should not af- 
fect severely the profiles of forbidden emission lines, so that 
they can be used to trace the geometry of the ejecta (but 
see Section H ] for the consequence of residual optical depth). 
IMaeda et al. | 2008j ) employed a more stringent criterion, 
restricting their sample to spectra with epochs 200 d to 
avoid any possible deformation of lines by optical-depth ef- 
fects. This criterion would reduce our sample from 98 to 53 
spectra (including 16 SNe not analysed by Maeda et al.), and 
deprive us of the possibility to investigate when the ejecta 
become fully transparent, which is addressed in Section [3] 

Applying a phase cut required a fairly precise estimate 
of the epoch of explosion or maximum light. Whenever a 
complete light curve was not available, this information was 
reconstructed from discovery and classification remarks re- 
ported in IAU-circulars. In exception to this rule, spectra 
of three SNe were included for which no e arly observations 
exist . However, the spectrum of SN 1995bb ()Matheson et al.1 
l200ll ) is decidedly nebular, as are the later two out of three 
spectra of SN 1990aj. The first one in this series appears 
rather peculiar and may still show some photospheric fea- 
tures, but was included for completeness. Finally, a spec- 
trum of SN 2005N was dated to ~ 90 d past maximum 
light by cross-correlation wit h a set of comparison spec- 
tra ()Harutvunvan et alJ 120081 ) . Besides the constraints on 
the phase, also spectra with insufficient signal-to-noise ratio 
(S/N) in the wavelength region of interest were rejected. 



2.1 The full sample 

The full catalogue of spectra studied in this work is pre- 
sented in Table [T] complemented by additional information 
on the SN classifications and host-galaxy properties. 



2.2 Previously unpublished spectra 

Our sample contains 26 nebular spectra of 17 SNe Ib/c 
not previously published elsewhere. Another 4 spec tra of 
the Asiago archive were shown by iTurattd (120031 ) and 
IValenti et al.l l|200Sl ) before. Most of these spectra were taken 
in the course of the ESQ- As iago SN monitoring programme 
in the 1990s (|Turattdl200ol ) using the ESO - La Silla 3.6m 
(equipped with EFOSC /EFOSC2), 2.2m (+ EFOSC2) and 
1.5m (+ Boiler & Chivens spectrograph) Telescopes and 
the 1.54m Danish Telescope (equipped with DFOSC). From 
2004 onwards, several spectra were acquired through ded- 
icated VLT programmes (VLT-U1 equipped with FORS2). 
The set is complemeted by single spectra taken with the Sid- 
ing Spring 2.3m Telescope (+ double-beam spectrograph) 
and the Nordic Optical Telescope (+ ALFOSC). Details on 
the dates of the observations and the instrumental setup are 
summarised in Table [2] 

Applying the selection criteria mentioned above, our 
full sample consists of 39 SNe with 98 nebular spec- 
tra. It contains almost all suitable spectra up to the 
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Table 1. List of SNe Ib/c included in the sample - 1 st part. 
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SN 


Type 


Host galaxy 


Morphology" 




Red shift 


Date 


Epoch c 


1983N 


lb 


NGC 5236 


SBc 


554 ±119 


0.0018(04) 


1984/03/01 


226 zb 5 


1985F 


Ib/c 


NGC 4618 


SBm 


544 ± 59 


0.0018(02) 


1985/03/19 


280 ± 4 


1987M 


Ic 


NGC 2715 


SABc 


12 16 zb 132 


0.0041(04) 


1988/02/09 


141 zb 7 














1988/02/25 


157 zb 7 


1988L 


Ic 


NGC 5480 


Sc 


1963 zb 100 


0.0065(03) 


1988/07/17 


90 zb 12 














1988/09/15 


149 zb 12 


1990B 


Ic 


NGC 4568 


Sbc(M) 


2255 zb 153 


0.0075(05) 


1990/04/19 


91 zb 2 














1990/04/30 


102 zb 2 


19901 


lb 


NGC 4650A 


SO /a(M) 


2880 zb 99 d 


0.0096(03) 


1990/07/26 


90 zb 2 














1990/12/21 


237 ± 2 














1991/02/20 


298 zb 2 


1990U 


Ic 


NGC 7479 


SBbc 


2525 zb 162 


0.0084(05) 


1990 /10 /20 


100 zb 12 
















104 zb 12 














1990/11/23 


134 zb 12 














1990/11/28 


139 zb 12 














1990/12/12 


153 zb 12 














1990/12/20 


161 zb 12 














1991/01/06 


178 zb 12 














1991/01/12 


184 zb 12 


1990W 


Ib/c 


NGC 6221 


SBc 


1481 ± 126 


0.0049(04) 


1991/02/21 


183 zb 3 














1991/04/21 


242 zb 3 


1990aa 


Ic 


MCG+05-03-016 


Sb 


5032 zb 108 


0.0168(04) 


1991/01/12 


130 zb 7 














1991/01/23 


141 zb 7 


1990aj 


Ib/Cp ec 


NGC 1640 


SBb(R) 


1604 ± QA d 


0.0053(02) 


1991/01/29 


140 ± 50 














1991/02/22 


164 zb 50 














1991/03/10 


180 zb 50 


1991A 


Ic 


IC 2973 


SBcd 


3232 ± 84 


0.0107(03) 


1991/03/22 


99 zb 10 














1991/04/07 


115 zb 10 














1991/04/16 


124 zb 10 














1991/06/08 


177 zb 10 


1991L 


Ib/c 


MCG+07-34-134 


Sc(M) 


9186 ± 200 d 


0.0306(07) 


1991/06 /08 


121 zb 20 


1991N 


Ic 


NGC 3310 


SABb(R) 


1071 ± 80 


0.0036(03) 


1991/12/14 


274 zb 15 














1992/01 /09 


300 zb 15 


1993J 


lib 


NGC 3031 


Sab 


— 140 zb 192 d 


-0.0001(06) 


1993/10/19 


205 zb 3 














1993 /11/19 


236 zb 3 














1993/12/08 


255 ± 3 














1994/01/17 


295 zb 3 














1994/01 /21 


299 zb 3 














1994/01 /22 


300 zb 3 














1994/03 /25 


362 zb 3 














1994/03 /30 


367 zb 3 


19941 


Ic 


NGC 5194 


Sbc(M) 


493 ± 70 


0.0016(02) 


1994/07/14 


97 zb 1 














1994/08 /04 


118 zb 1 














1994/09 /02 


147 zb 1 


1995bb 


lb /c 




S /Irr 


1626 zb 250 


0.0054(08) 


1995/12/17 




1996D 


Ic 


NGC 1614 


SBc(M) 


4531 zb 167 


0.0151(06) 


1996 /09/10 


214 zb 10 


1996N 


lb 


NGC 1398 


SBab(R) 


1396 zb 214 


0.0047(07) 


1996/10/19 


224 zb 7 














1996 /12/16 


282 zb 7 














1997/01/13 


310 zb 7 














1997/02/12 


340 zb 7 


1996aq 


lb 


NGC 5584 


SABc 


1675 zb 83 


0.0056(03) 


1997/02/11 


176 zb 4 














1997 /04/02 


226 zb 4 














1997/05/14 


268 zb 4 


1997B 


Ic 


IC 438 


SABc(R) 


2919 zb 144 


0.0097(05) 


1997/09/23 


262 zb 5 














1997/10/11 


280 ± 5 














1998/02/02 


394 zb 5 


1997X 


Ic 


NGC 4691 


SBO/a 


1072 zb 47 


0.0036(02) 


1997/05/10 


103 zb 5 














1997/05/13 


106 zb 5 


1997dq 


Ic 


NGC 3810 


Sc 


993 zb 114 d 


0.0033(04) 


1998/05/30 


217 zb 10 














1998/06/18 


236 =b 10 


1997ef 


BL-Ic 


UGC 4107 


Sc 


3452 =b 63 


0.0115(02) 


1998/09/21 


287 =b 3 


1998bw 


BL-Ic 


ES0184-G82 


SBbc 


2445 zb 135 


0.0082(05) 


1998/09/12 


126 zb 1 














1998/11/26 


201 zb 1 














1999/04/12 


337 =b 1 














1999/05/21 


376 zb 1 


1999cn 


Ic 


MCG+02-38-043 


S a l:i 


6502 zb 235 


0.0217(08) 


2000/04/08 


297 =b 5 


1999dn 


lb 


NGC 7714 


SBb(M) 


2744 zb 74 


0.0091(02) 


2000/09/01 


375 zb 5 



Gaskcll et al. 1986 
Filippcnko & Sargent 1986 
FilippenkoetaL 1990 
^iligp^nkoe^L^990 
MathesonetaL 2001 
^atheson^^a^^OOl 
^tocchiatt^e^L^OOl 
MathesonetaL 2001 
JSlmhaindi^e^^^OM 
ElmliamdietaL 2004 
BlmliamdietaL 2004 
^athesrajj^a^^OOl 
^^n^^^^^^p^^^^94 
Asiago archive 
Matheson ct al. 2001 
MathesonetaL 2001 
Asiago archive 
^atheson^^a^^OOl 
Gomez &:^L6p ^ez^^l994 



Asiago archive 
Asiago archive 
Gomez fc_L££ ez_2002 



Matheson ct al. 2001 

Asiago archive 

Asiago archive 

MathesonetaL 2001 

Gomez &c Lopez 1994 
^atheson^^a^^OOl 
^^n^z^^^L^^z < ^^94 

Gomez & Lopez 1994 

Gorasz&jjojez 2002 
^atheson^^a^^OOl 

Matheson et al. 2001 



Barbon ct 


al. 


1995 


Barbon ct 


al. 


1995 


Barbon ct 


al. 


1995 


Barbon ct 


al. 


1995 


Barbon ct 


al. 


1995 


Barbon ct 


al. 


1995 


Barbon ct 


al. 


1995 


Barbon ct 


al. 


1995 



Filippcnko ct al. 1995 
^lip£enko^^L^95 
Filippenkoe^L^95 
Mathesone^h 2001 
Asiago archive 
Sollcrman ct al. 1998 
Sollcrman ct al. 1998 
SoUermane^L 1998 
Sollermaiie^L 1998 
Asiago archive 
Asiago archive 
Asiago archive 
Asiago archive 
Asiago archive 
Asiago archive 
Gome^Uj6gez^002 
Asiago archive 
Asiago archive 
Matheson et al. 2001 



Mat heson et al. 2001 



Patat 


ct 


al. 


2001 


Patat 


ct 


al. 


2001 


Patat 


ct 


al. 


2001 


Patat 


ct 


al. 


2001 



Asiago archive 
Asiago archive 



year 2004 that could be retrieved from the litera- 
ture, complemented by previously unpublished spectra 
from the Asiago Supernova Archive (jBarbon et all 1 19931 . 
http://web.oapd.inaf.it/supern/cat/) and selected spectra 
obtained through dedicated programmes after 2004. The ob- 
servations thus span the entire era of SN CCD spectroscopy. 



Spectra from the Asiago and MPA archives a re pre- 
sente d in Fig. [1] They have been optimally extracted ijHornel 
1986) using standard tasks in iRArQ or MIDAS, wavelength 



1 I RAF is distributed by the National Optical Astronomy Obser- 
vatories, which are operated by the Association of Universities 
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Table 1. continued. List of SNe Ib/c included in the sample - 2 nd part. 



SN 


Type 


Host galaxy 


Morphology a 


v* eo 


Redshift 


Date 


Epoch c 


Reference 


2000cw 


Ic 


NGC 3810 


Sc 


1049 zb 114 


0.0035(04) 


2001/03/17 


112 zb 12 


A. si ago archive 


2002ap 


BL-Ic 


1\ Va \^/ vjZO 


Sc 


657 zb 22^ 


n 0099/01 

U . U vl z z \ u ± j 


9009 /Of-i /OQ 

ZVfVIZ [ VIV3 f VI .3 


123 zb 1 


Foley et al. 2003 














9009 /(\(\ /is 

Z, U U A j U U / to 


132 zb 1 


PTnloir of al 900 1 ? 














9009 /07 / 1 9 

ZVIVf Z / VI { 1 ^-Z 


156 zb 1 


T7nl 1317 tat al 900^ 














9009 /OS /OQ 




Jroicy et ai. zvivio 














9009 /I O /01 


237 zb 1 


roicv et ai. zuua 














9009 /1 /OQ 


245 zb 1 


Foley et al. 2003 














2002/10 /14 


250 zb 1 
















2002/11/06 


274 ± 1 


Foley et al. 2003 














2003/01/07 


336 zb 1 


Foley et al. 2003 














2003/02/27 


386 zb 1 


Foley et al. 2003 


2003jd 


BL-Ic 


MCG-01-59-021 


SABm 


5654 zb 78 


0.0188(03) 


2004/09/11 


317 zb 1 


Mazzali et al. 2005 














900zi /in/is 


354 zb 1 


a/To 1710 1 ; o+ ^^ 90o>^ 


2004aw 


Ic 


NGC 3997 


SBb(M) 


4900 zb 118 


0.0163(04) 


2004/11/14 


236 zb 1 
















2004/12/08 


260 zb 1 


Taub enb crger et al. 20 C) 6 














2005/05/11 


413 zb 1 


Taubenberger et al. 2006 


2004gt 


Ic 


NGC 4038 


SBm(M) 


1424 ± 133 


0.0047(04) 


2005 /05 /24 


160 zb 5 


Asiago archive 


2005 N 


lb /c 


NGC 5420 


Sb 


4885 ± 198 d 


0.0163(07) 


2005 /01/21 


88 zb 30 


H ar ut y uiiy an ct al . 2008 


_)()| ( (j p 


lb 


in \j yo J 


Sc(M) 


4- 1 sn 
4: / \j zn iou 


vi ,\j ±-±z\vivt j 


900fi /1 1 /1 fi 

ZVIVJVJ 1 -L-L/ 1U 


312 zb 7 


Macda ct al 2008 


2006T 


lib 


NGC 3054 


SBb(R) 


2560 zb 182 


0.0085(06) 


2007/02/18 


371 zb 2 


Macda ct al. 2008 


2006aj 


BL-Ic 






9845 zb 250 


0.0328(08) 


2006/09/19 


204 zb 1 


Mazzali et al. 2007a 














2006/11/27 


273 zb 1 


MPA data base 














2006/12/19 


295 zb 1 


MPA data base 


2006gi 


Ic 


NGC 3147 


Sbc 


2820 zb 178 d 


0.0094(06) 


2007/02/10 


148 ± 5 


Asiago archive 


20061d 


lb 


UGC 348 


SABd(R) 


4168 zb 41 


0.0139(01) 


2007/07/17 


280 ± 4 


MPA data base 














2007/08/06 


300 ± 4 


MPA data base 














2007/08/20 


314 =b 4 


MPA data base 


2007C 


lb 


NGC 4981 


SBbc(R) 


1766 zb 116 


0.0059(04) 


2007/05/17 


131 ± 4 


MPA data base 














2007/06/20 


165 =t 4 


MPA data base 


20071 


BL-Ic 


anonymous 


late spiral 


6445 zb 250 


0.0215(08) 


2007/06/18 


165 =b 6 


MPA data base 














2007/07/15 


192 ± 6 


MPA data base 



a Classification according to LEDA (Lyon-Mcudon Extragalactic Database, http://lcda.univ-lyonl.fr/). 

b Recession velocity in kms" 1 , inferred from narrow Ha: error from ' vmaxg ' (LEDAJ. 

c Epoch in days from S-band maximum light; photometry of SN 1985F from lTsvetkovl Jl986fi . 

No Hck visible; heliocentric host-galaxy recession velocity from NED (NASA/IPAC Extragalactic Database, http://nedwww.ipac.caltech.edu/) used. 



calibrated with respect to arc lamps, and flux calibrated us- 
ing instrumental response curves obtained from spectropho- 
tometric standard stars observed in the same nights. How- 
ever, no attempt has been made to calibrate the fluxes to a 
proper absolute scale through a comparison to contempora- 
neous photometry. 

The spectra all show [O i] AA6300, 6364 as one of 
their strongest features, complemented by other lines typ- 
ical of SNe Ib/c at late phases, most notably [Call] 
AA7291,7323/[On] AA7320, 7330, Mgi] A4571, a multitude 
of blended [Fell] lines around 5000 A, and the near-IR Call 
triplet. In spectra taken less than 150 d after maximum a 
contribution from photospheric lines and a weak pseudo- 
continuum can be discerned. Some spectra also show evi- 
dence of an underlying stellar continuum caused by an im- 
perfect host-galaxy subtraction. 



3 FITTING THE OXYGEN LINE 

[O i] AA6300, 6364 is one of the dominant features of nebular 
spectra of stripped-envelope CC-SNe, and the most useful 
to study the ejecta geometry, in particular the degree of as- 
phericity. Compared to the multitude of forbidden Fe lines 
found mostly at bluer wavelength and to the [Ca n] / [O n] 
feature around 7300 A, the [Oi] AA6300, 6364 doublet is 

for Research in Astronomy, Inc, under contract with the National 
Science Foundation. 



largely isolated and unblended. In contrast to the weaker 
Mg i] A4571 it lies in a region which is covered by almost all 
late-time spectra, and where the sensitivity of most spectro- 
graphs is at its maximum, allowing for relatively good S/N. 
Furthermore, oxygen is the most abundant element in the 
ejecta of stripped-envelope CC-SNe, thus tellingmore about 
the overall geometry than the distribution of a minor species 
like Ca. 

In this section we describe our method to obtain con- 
straints on the geometry of the ejected oxygen, which in- 
volves the characterisation of the expected line profile, the 
development of a suitable parametrisation, and the actual 
procedure applied to infer the best-fitting parameters for 
each spectrum. 

3.1 Line profiles 

In completely transparent SN ejecta the profile of a forbid- 
den emission line traces the emissivity in this line, which in 
turn is determined by the spatial distributions and velocity 
fields of both the emitting species and 56 Co, whose radioac- 
tive decay provides the energy to e xcite the line's upper level 
l|Fransson fc Chevalier|[l987lTl989l '). In SNe Ic the O- and Co- 
rich parts represent a significant fraction of the entire ejecta, 
and the [O i] feature traces a substantial amount of material. 
Thanks to homologous expansion (r = vt) , the profile of the 
[Oi] emission line is a ID line-of-sight projection of the 3D 
oxygen emissivity distribution. 

In this work, we fit the oxygen feature with a Gaus- 
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Figure 1. Spectra from the Asiago archive and the MPA data base, most of them previously unpublished (cf. Table |3J. All spectra are 
shown at their rest wavelength inferred (whenever possible) from narrow interstellar Ha lines. Overly strong continuum slopes and narrow 
host-galaxy emission features have been removed for presentation purposes, and the spectra have been scaled and vertically displaced by 
arbitrary amounts. The major features are labelled in the first spectrum. 



sian or - for SNe with more complex line profiles - with a 
superposition of multiple Gaussians. This provides at least 
qualitative information on the distribution (through the po- 
sition and strength of the various components) and radial 
extent (through the components' FWHM) of excited oxy- 
gen in the SN ejecta. However, a full restoration of the 3D 
density distribution is not attempted, since the solution is 
highly degenerate. 



3.2 [OI] A6300 and [OI] A6364 

As mentioned above, an advantage of the oxygen feature 
is its isolated position, unblended with lines from other 
elements. However, the feature itself is a doublet of [O i] 
A6300 and [O i] A6364, both forbidden Ml transitions which 
share the same upper level ( 3 Pi,2 - 1 ^>2)- The intensity 
ratio of these two lines depends on the ambient O I den- 
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Table 2. Instrumental details of spectra from the Asiago archive 
and the MPA data base. 



SN 


Date 


Instrumental setup 






1990U 


1990/11/23 


ES0 3.6m 


+ EFOSC 


f B300 + 


R300 




1990/12/20 


ES0 3.6m 


+ EFOSC 


+- B300 




1990W 


1991/02/21 


ES0 3.6m 


+ EFOSC 


f B300 + 


R300 




1991/04/21 


ES0 3.6m 


+ EFOSC ■ 


+■ B300 + 


R300 


1990aj 


1991/01/29" 


ES0 2.2m 


+ EFOSC2 


+ grl 






1991/02/22" 


ES0 3.6m 


+ EFOSC 


f B300 + 


R300 


1996D 


1996/09/10 


ESO 1.5m 


+ B&C + gtl5 




1996aq 


1997/02/11 


ES0 3.6m 


+ EFOSC ■ 


f R300 






1997/04/02 b 


ESO 1.5m 


+ B&C + i 


»tl5 






1997/05/14 


ESO 2.2m 


+ EFOSC2 


+ grl + 1 




1997B 


1997/09/23 


ESO 2.2m 


+ EFOSC2 


+ gm5 






1997/10/11 


Danish 1.54m + DPOSC + gr5 






1998/02/02" 


ESO 3.6m 


+ EFOSC2 


+ gr6 




1997X 


1997/05/13 


ESO 2.2m 


+ EFOSC2 


+ grl + 1 


Sr5 


1997dq 


1998/05/30 


ESO 3.6m 


+ EFOSC2 


+ B300N 


+ R300N 


1999cn 


2000/04/08 


ESO 3.6m 


+ EFOSC2 


+ grll 




1999dn 


2000/09/01 


ESO 3.6m 


+ EFOSC2 


+ grl2 




2000ow 


2001/03/17 


Danish 1.54m + DPOSC + gm4 




2002ap 


2002/10/14 


SSO 2.3m ■ 


f DBS 






2004gt 


2005/05/24 


VLT-U1 + 


FORS2 + 


300V 




2006aj 


2006/11/27 


VLT-U1 + 


FORS2 + 


300V + 3001 




2006/12/19 


VLT-U1 + 


FORS2 + 


300V 




2006gi 


2007/02/10 


NOT 2.56m + ALFOSC + gm4 




20061d 


2007/07/17 


VLT-U1 + 


FORS2 + 


300V 






2007/08/06 


VLT-U1 + 


FORS2 + 


300V 






2007/08/20 


VLT-U1 + 


FORS2 + 


300V 




2007C 


2007/05/17 


VLT-U1 + 


FORS2 + 


300V 






2007/06/20 


VLT-U1 + 


FORS2 + 


300V + 3001 


20071 


2007/06/18 


VLT-U1 + 


FORS2 + 


300V + 3001 




2007/07/15 


VLT-U1 + 


FORS2 + 


300V 





" Already shown bv lTurattol )2003t). 

6 Already shown bv lValenti et al.l j20ol) . 



sity, and can vary from 1:1 to 3:1 depending on the en- 
vironmental conditions. The transition between the asymp- 
totic values occurs at Ol densities of n(Ol) ~ 10 10 cm" 3 . 
This fact has been theo retically derived by iLi fc McCra^l 
lll992h andlChueail (Il992ri and obscrv ationally confirmed by 
iPhillips fe Williams! (Il99ll) and ISpvromilio fc Pintol (|l99lh 
for SN 1987A, iLeibundgut et al l (|l99lh for SN 1986J. and 
ISpvromiliol l|l99ll ) for SN 1988A. 

Compared to the aforementioned SNe II, the stripped 
CC-SNe of our sam p le sho w larger ejecta velocities. While 
ISpvromilio fc Pintol (|l99ll ) measured a FWHM of 2800 
km s" 1 for the [Oi] line in SN 1987A, the SNe discussed 
here have FWHM of ~ 6000 kms -1 . Assuming a Gaussian 
density profile with 6000 kms"" 1 FWHM and 1M Q of neu- 
tral oxygen homogeneously distributed within the ejecta, the 
central O I density would have dropped to 4-5 x 10 8 cm" 3 
by 100 d, an order of magnitude below the density where de- 
viations from a 3 : 1 line ratio become apparent. Even if the 
oxygen was clumped on small scales wit hin an overall Gaus- 
sian p rofile [as in SN 1987A, for which ISpvromilio fc Pintol 
l|l99lf ) suggested an oxygen filling factor of ~ 10% based on 
the observed evolution of the [O i] ratio] , ratios significantly 
different from 3 : 1 would only be encountered for very small 
filling factors (< 10%). We therefore adopt the low-density 
limit (3 : 1) for all spectra, noting that possible small devia- 
tions at the earliest epochs do not severely affect any of our 
basic conclusions. 

With this choice, it is sufficient to specify the ampli- 
tude, central wavelength and FWHM of the A6300 line. The 
parameters for the A6364 line are then forced. Hereafter we 



refer to such a set of three parameters as one component of 
the fit to the [O i] profile; up to three such components were 
employed to obtain good fits. 

3.3 The fitting procedure 

The actual fitting was accomplished using the iraf task 
NFITlD, which is part of the STSDAS package. A user-defined 
fitting-function was introduced, which consisted of up to 
three components, as defined above. 

The background level was determined by eye on both 
sides of the [O i] feature and interpolated linearly to ac- 
count for a possible underlying continuum formed by resid- 
ual photospheric lines or a contamination by star light. Since 
the background parameters are determined from a different 
wavelength region than the line parameters, the background 
fit is technically decoupled from the line fit, although of 
course the best-fit line parameters may be affected by the 
choice of the background. This can be problematic in rela- 
tively early spectra (up to ~150d), where underlying pho- 
tospheric lines are present, and the local continuum is more 
strongly inclined than at later times. However, tests with 
different choices of the continuum level have shown that the 
uncertainty introduced e.g. in the central wavelength of the 
line is <J 5 A even in cases with very complex background. 
More typical uncertainties are < 2 A and thus smaller than 
the uncertainties in the wavelength calibration and redshift 
correction of most spectra. 

For the up-to-nine parameters of the line fit (amplitude, 
position A and FWHM of up-to-three components) the main 
difficulty consisted in identifying the global minimum for a 
given number of components. Therefore, for fits with two or 
three components, a refined two-step fitting procedure was 
applied to avoid local minima. First, a set of synthetic line 
profiles were generated, changing all parameters by equidis- 
tant steps over a fairly wide range. The resulting profiles 
were compared to the observed ones, and the one with min- 
imum RMS was identified. In a second step the values thus 
derived were used as initial guesses for nfitID, ensuring 
that the fit converged to the global minimum. We always 
started the fitting with one component. Other components 
were added only if the fit residuals strongly exceeded the 
noise level of the spectrum. 

The best-fitting values for single- and double- 
component fits to all spectra are reported in Appendix 1X1 
( Table lAl) . while the parameters of three-component fits to 
individual spectra are given in Table lA~2l In these tables and 
during the further discussion, we define as the integrated 
flux of the i-th component of the fit, normalised to the to- 
tal integrated flux of the oxygen feature. The best fits are 
compared with the observed line profiles in Fig. [2] 



4 BLUESHIFTED LINE CENTROIDS AT 
EARLY PHASES 

In Fig. [3] we plot the position of the A6300-Gaussian in the 
one-component fits against the epoch of the spectra. The 
position should be a fair tracer of the actual line centroid. 
The scatter of the data points arises from the peculiarities 
of individual objects and uncertainties in their redshifts. 
However, on top of this scatter Fig. [3] shows that there 
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Figure 2. Observed [Ol] AA6300, 6364 features in the spectra of our sample (solid black lines) with multi-Gaussian fits (Tables lAll and 
I A2t overplotted (solid red lines). Individual components are overplotted as dotted blue lines, and the adopted linear background levels 
are indicated by dashed grey lines. Ha was subtracted from the spectrum of SN 1993 J (cf. Table lAl) . The different line categories as 
defined in Scction l5.1l are labelled as follows: GS for Gaussian, NC for narrow core, DP for double peak, AS for asymmetric / multi-peaked 
(alternative classifications are given in brackets). Only one spectrum of each object is included in the figure. 
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150 200 250 300 350 
days from B bond maximum 

Figure 3. [O i] A6300 line centroids (as inferred from the position 
of the A6300 Gaussian in one-component fits), plotted vs. the 
epochs of the spectra with respect to B-band maximum. The 
filled circles represent bins of 10 spectra. A systematic blueshift 
can be discerned at epochs earlier than 200 d. 



identified components 
per wavelength bin; 
weight =1 for each SN 



I I all spectra 
M<t only <1 50d 




6220 6240 6260 6280 6300 6320 6340 6360 6380 
wavelength [A] 

Figure 4. Histogram of components identified in the spectra with 
our Gaussian fitting procedure (wavelengths refer to [O i] A6300; 
cf. Tables Ell E2]and Fig.E}. To give equal weight to all SNc, the 
numbers have been rescaled such that every SN yields a contribu- 
tion equivalent to one component (hence the fractional numbers). 
The empty histogram refers to the full sample, the shaded region 
to a subsample of spectra taken earlier than 150 d after maximum. 



is a systematic trend of the [O i] feature being blueshifted in 
spectra taken earlier than ~ 200 d past maximum, and the 
effect is stronger the earlier the phase. In the spectra taken 
around 100 d, the average blueshift is ~20A, corresponding 
to ~ 1000 kms -1 . In the following, possible interpretations 
of the observed effect are discussed, and their suitability to 
explain the observations is considered. 

(i) Ejecta geometry. Line shifts such as those observed 
in [O i] could, in principle, arise from a one-sided ejecta ge- 
omet ry, caused for instance by l ow-mode convective i nstabil- 
ities (jScheck et al ] |2004l . 120061 : iKifonidis et ai][ 2006) or the 
Stan ding Accretion Shock Instability (SASI; iBlondin et al.l 
2003). However, this can not explain the decrease of the 
shifts with time. Moreover, there is no reason why particu- 
lar ejecta geometries should result in a systematic 6/tteshift 
of the line centroid, as different spatial orientations should 
occur in the observed sample. Note that relativistic forward 
boosting is irrelevant at the observed ejecta velocities (no 
more than ~8000 kms -1 , even in the extreme line wings). 

(ii) Dust formation. As the SN ejecta expand and cool, 
the temperature eventually drops below the threshold where 
dust can form. Consequently, the light from the far side 
of the ejecta is partly absorbed, resulting in the suppres- 
sion of the redshifted part of emission lines. This effect has 
been observed in some Type II SNe, for instance SNe 1987A 
and 1999em, more than one year aft er e xplosion (cf. e.g. 
Danzi ger et all 1 1989 . iLucv et al ] ll989l and LEhnhamdi et al.l 
2003). In ordinary SNe Ib/c dust formation has never 
been unambiguously detected at a few hundred days (e.g.. 



Sollerman. Leibundeut fc Spvromiliolll998l . 
2000 and Elmhamdi et al .1120041 ; but see also 



Mathcson ct al 



Matthews et al 



20021 ). Moreover, dust formation (if present) should manifest 
itself in a line blueshift increasing with time, the opposite of 
what we see in our sample, and hence cannot be a suitable 
explanation for our observations. 

(iii) Contamination from other emission lines. In prin- 
ciple, other lines blended into the blue wing of [O i] 
AA6300, 6364 could generate the observed blueshift of this 



feature. At early epochs, when a strong blueshift is ob- 
served, the contamination could e.g. arise from residual 
emission of permitted photospheric lines. The biggest prob- 
lem w ith this idea is t h e app arent lack of suitable candi- 
dates. lElmhamdi et all l|2004l ) suggested a contribution of 
Fen A6239 in SN 19901 around +90 d, but this is not ex- 
pected to be a particularly strong line, and it is unclear why 
it should be so prominent while other, intrinsically stronger 
Fe lines are not. Moreover, a histogram of fit components 
(Fig. |4]| does not show evidence of a distinct additional line 
at a wavelength shorter than 6300 A. Instead, in the sub- 
sample of spectra taken at < 150 d (shaded area in Fig. [4]), 
the distribution of fit components smoothly smears out to 
shorter wavelengths. Finally, in Section [7] we will show that 
the profile of Mgi] A4571 is similar to that of [O i] A6300 in 
a majority of our spectra (also those with blueshifted lines), 
and an identical contamination in both lines is very unlikely. 

(iv) Opaque inner ejecta. The failure of other explana- 
tions and the characteristics of the fit-component histogram 
leave us with residual opacity in the core of the e jecta as 
a po s sible explanation for the observed blueshift l|Chugail 
ll992l : IWang fc Hulll994h . Optically thick inner ejecta could 
prevent light from the rear side of the SN from penetrating, 
creating a flux deficit in the redshifted part of emission lines. 
The opacity could be caused by e.g. densely packed weak Fe 
transitions (electron-scattering turns out to be at least an 
order of magnitude too weak). To see the effect on the line 
profile, we created a simple model using a Monte Carlo code 
(see Fig. O, where photons are absorbed or scattered on 
their way to an observer with a probability proportional to 
the ambient matter density (grey opacity). A profile calcu- 
lated for an unrealitically early epoch of 30 d assuming pure 
elastic electron-scattering shows a characteristic tail on the 
red side, but too little blueshift of the line core to be con- 
sistent with observations at 100 d (see Fig. [SJ). This demon- 
strates that Thomson scattering provides too little opacity, 
and, furthermore, would modify the line profile in an un- 
desired way if it were strong enough. If, instead, the calcu- 
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Figure 5. Synthetic profiles of [Ol] AA6300, 6364 calculated us- 
ing a Monte Carlo code. A Gaussian density distribution, and 
both an emissivity and opacity proportional to the density have 
been assumed. The pure-scattering calculation is based on the 
Thomson cross-section crxh for e — -scattering, assuming on aver- 
age singly ionised material. The cases of pure absorption were 
computed using a grey opacity with a = 10<t t1i . For comparison, 
the unabsorbed profile is also shown. 



lations are performed for grey absorption, good results can 
be obtained if the cross section is chosen appropriately. In 
particular, the observed time evolution of the line blueshift 
is reproduced qualitatively thanks to the t~ 2 scaling of the 
column density. 



(ii) Narrow-core lines (NC), characterised by an addi- 
tional narrow component centred close to the rest wave- 
length and atop the broad base of the line (offsets <C 20 A). 
These are quite frequent and can be explained (a) by ax- 
isymmetric explosions with the emitting oxygen located in 
a t orus or disk perpendicular to the line of sig ht (as inferred 
bvlMazzali et aUldooj . iMaeda et al.ll2002l and lMaeda et all 
2006 for SN 1998bw), (b) by spherically symmetric ejecta 
with an enhanced core density, or (c) by a blob of oxygen 
moving nearly perpendicularly to the line of sight [cf. (iv)]. 
Note that the presence of dense cores has been suggested 
for several SNe to explain their line profiles and late-time 
light-curve slopes lllwamoto et al. 2000l; Mazzali et al.ll20o"ol ; 



IMaeda et al.ll2003l ; IMazzali et aLlkoOTal lbfT 



(iii) Double-peaked profiles (DP) with two compara- 
bly strong components, one blueshifted and the other red- 
shifted by similar amounts. These are most readily ex- 
plained by a torus-shaped oxygen distribution viewed nearly 
sideways (from ang les of ~ 60-90° to th e symmetry axis; 
IMazzali et~aDl2005l and lMaeda et alJliooi ). No double peak 
can be realised in spherical symmetry. Hence, this class of 
line profile re quires asphericity. The prototype of thi s class 
is SN 2003jd ijMazzali et al.ll2005l : IValenti et al.ll200ct ). 

(iv) Multi-peaked or asymmetric profiles (AS), pro- 
duced by additional components of arbitrary width and 
shift with respect to the main component. These profiles 
are either indicative of ejecta with large-scale clumping, a 
single massive blob, or a unipolar jet. Like the double peaks, 
they cannot be reproduced within spherical symmetry. 



5 STATISTICAL ANALYSIS, INFERRED 
EJECTA GEOMETRIES 

The main intention of the multi-Gaussian line fitting is to 
derive information on the occurrence of different ejecta ge- 
ometries in the sample. Of course, without additional as- 
sumptions it is not possible to restore the full 3D density 
distribution from its ID projection given by the line profiles. 
While a forward calculation of emerging profiles for a given 
density distribution is staightforward, backward inference is 
highly degenerate. 

A brief overview of some possible ejecta geometries and 
the corresponding observed line profiles is given in Table |3J 
illustrating that in several cases the geometry cannot be 
determined with confidence. Yet, for most observed profiles 
we can at least exclude certain configurations. 

5.1 Taxonomy 

Guided by the results of the multi-Gaussian fitting, we in- 
troduce four principal classes of line profiles. Note that this 
classification scheme is a simplistic choice, based on experi- 
ence acquired during the fitting. 

(i) Gaussian profiles ( GS), well reproduced by single- 
component fits, with the residuals showing no evidence of a 
second component within the noise level. These profiles are 
expected from spherically symmetric ejecta with a nearly 
Gaussian emissivity distribution, but could alternatively be 
the outcome of e.g. axisymmetric explosions viewed from 
intermedi ate angles (40-50°, depending on the degree of as- 
phericity; IMaeda et aT]|2008r ) . 



Assigning our sample of SNe to these categories is some- 
times ambiguous. For instance, an [Ol] feature which con- 
sists of a main peak and a Doppler-shifted blob may appear 
double-horned, and can be confused with genuine double 
peaks formed by a toroidal oxygen distribution as defined 
in (iii). A further complication for the classification arises 
from possible bulk shifts of the [O i] feature with respect to 
its rest wavelength: as discussed in Section \$\ this does not 
necessarily have a geometric origin. 

Fig. [^indicates the class membership of individual SNe. 
Only SNe fitted with two components [classes (ii) to (iv)] 
are included in this figure, which shows the absolute wave- 
length offset between the two fit components (Table lATjl as 
a function of a w , the normalised flux of the weaker com- 
ponent. In this diagram, narrow-core SNe [class (ii)] popu- 
late a strip along the abscissa (|Ai — Aa | < 20 A), double 
peaks [class (iii)] an area of larger wavelength offset and 
0.4 <C q w 0.5, while SNe with asymmetric or multi-peaked 
profiles [class (iv)] are mostly contained in a region charac- 
terised by a w <J 0.3 and |Ai — A2 > 20 A. Note, however, 
that some objects lying in the narrow-core strip actually 
belong to class (iv), since their weaker component has too 
large an offset (> 20 A) from the rest wavelength to fulfil 
the criteria defined for class (ii). 



5.2 Statistical evaluation 

In Table [3] the statistical summary of this analysis is pre- 
sented. As we will see below, deviations from spherical sym- 
metry affect all types of stripped-envelope CC-SNe, and 
are not reserved to particularly energetic or highly-stripped 
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Table 3. Selected oxygen geometries and corresponding line profiles. 



Oxygen emissivity distribution Line profile Global symmetry 



radial Gaussian 


Gaussian 


spherically symmetric 


enhanced central density 


narrow core on top 


spherically symmetric 


hard-edged homogeneous sphere 


parabolic 


spherically symmetric 


spherical shell 


flat-topped 


spherically symmetric 


torus viewed from top 


narrow core 


axisymmetric 


torus viewed from the side 


double peak, symmetric to Ao 


axisymmetric 


torus viewed from intermediate angle 


Gaussian-like 


axisymmetric 


small-scale dumpiness 


fine-structured peak 


asymmetric 


unipolar jet, one-sided blob 


extra-peaks / shoulders, 


asymmetric 




shifted with respect to Aq 
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Figure 6. Absolute wavelength difference between two components found by multi-Gaussian fitting of [Oi] A6300, as a function of a w , 
the relative flux of the weaker component. Objects with single-Gaussian line profiles (cf. Fig. are not shown. Filled blue symbols stand 
for narrow-core SNe, filled red symbols for double peaks, and open symbols for SNe with asymmetric or clumpy ejecta. The different 
classes appear fairly well separated in this plot. 



events . Thi s is in agreem e nt wit h the results of lModiaz et alj 
|200St ) and lMaeda et all (|200Sl ). 

Spherically symmetric objects. SNe whose [O i] profiles 
are well fit with single Gaussians make up little more than a 
quarter of all objects, even within the uncertainties. Consid- 
ering all possibly spherical SNe [i.e., classes (i) and (ii)], and 
again including objects with ambiguous classification, we 



find their fraction to be just over 50%. Given that for some 
of these objects the S/N is too low to identify more than 
one component, that also jet- like explosions yield single- 
peaked symmetric profiles if viewed not too far from the 
jet axis, and that blobs moving roughly perpendicular to 
the line of sight can mimic narrow line cores, it is evident 
that this is really an upper limit for the number of spher- 
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Table 4. SN taxonomy in terms of [Ol] line profiles. The er- 
rors account for possible alternative classifications as indicated in 
Fig. El 



Category Number Percentage 



(i) Gaussian 


7 +4 
1 -2 


18 


(ii) Narrow core 


"tl 


28 


(iii) Double peak 


6±i 


15 


(iv) Asymmetric / blobs 


15±S 


39 



ical obje cts. This is e xemp lified by SN 1997X [class (i)], 
for which I Wang et all lj200lJ l measured exceptionally strong 
continuum polarisation, of the order of 4%, clearly indicat- 
ing globally aspherical ejecta. However, at the same time 
the [O i] profile reveals no obvious deviation from spherical 
symmetry (Fig. [5]). Hence, SN 1997X is probably intrinically 
aspherical, but viewed from a direction in which the line-of- 
sight projection of the oxygen emissivity mimics a spherical 
explosion. We thus speculate that probably more than half 
of all stripped-envelope CC-SNe are significantly aspherical. 

Symmetric double peaks. Between 5 and 18% of the SNe 
in our sample belong to class (iii), i.e. their line profiles 
are best reproduced by a symmetric double-peak configu- 
ration. This suggests a somewhat lower occurrence ra te of 
double peaks t han t he samples of iMaeda et al.l (|2008r i and 
iModiaz et aD (|2008t l. where 28% and 37% of the SNe, re- 
spectively, were double-peaked. However, it sh ould be noted 
that some of the doub l e-peak ed objects of IMaeda et al.l 
l|2008h and lModiaz et alJ l|2008l h such as SNe 2004ao, 2005aj 
and 2005bf, seem to lack symmetry about Ao, and might 
have been placed in class (iv) in our scheme. 

Jet-SNe. In the jet-models of lMaeda et"aH (j2006l , l2008r i. 
oxygen is distributed in a torus-like geometry perpendicular 
to the jet axis, and the [O i] profile is strongly viewing-angle 
dependent (see Table [3}. These models predict the ratio of 
narrow cores to double peaks to be ~ 1 : 5, rather insensi- 
tive to the degree of asphericity. On the contrary, we find 
fewer double-peaked profiles than narrow cores, suggesting 
that the majority of these narrow cores do not originate 
from jets, but e.g. from an enhanced central density. This 
also means that only a rather small fraction of all SNe Ib/c 
have a bipolar-jet ge ometry. In fact, de pending on whether 
moderately (BP2 of IMaeda et all 120061 ) or strongly (BP8) 
aspherical models are used to evaluate the viewing-angle de- 
pendence of the profiles in detail, we find jet-SN fractions of 
~ 50% or ~ 25% in our sample. 

CC-SN subtypes. In Fig. [7]we show how the traditional 
stripped-envelope CC-SN subtypes (i.e., BL-Ic, Ic, lb and 
lib) are distributed in terms of line-profile classes. Although 
the total number of objects is too small for robust state- 
ments, we are tempted to attribute some significance to the 
trends we can discern. While SNe Ic, which form the ma- 
jority of our sample (59%), are relatively homogeneously 
distributed, SNe lb (21% of our sample) belong mainly to 
the multi-peaked / aspherical category. Finding objects with 
extended envelopes to show particularly strong aspericity 
appears counter-intuitive, and we have no convincing expla- 
nation for this behaviour. Similarly, broad-line SNe Ic (BL- 
Ic; ~ 15% of our sample) show a tendency towards nebular 




Figure 7. Allocation of the SNe in our sample to the line-profile 
classes defined in Section l5.ll GS stands for 'Gaussian' [class (i)], 
NC for 'narrow core' [class (ii)], DP for 'double-peaked' [class (iii)] 
and AS for 'asymmetric/ multi-peaked' [class (iv)] (cf. Table 3}. 
SNe lib, SNe lb, SNe Ic and BL-Ic SNe are separately shown in 
the different panels. 

[O i] lines with narrow cores. At first, this seems to support 
the jet model. However, as discussed above, if all broad-line 
SNe Ic with narrow core were interpreted as jet-events, a 
much larger number of double peaks would be expected|j 
The FWHM of [O i] A6300 (taken from the one-component 
fits), averaged over all spectra, increases from SNe Ib/IIb 
(5205 ± 862 kms" 1 ) over normal-energetic SNe Ic (5942 ± 
1376 kms" 1 ) to broad-line SNe Ic (7343 ± 1724 kms" 1 ), 
reflecting the trend found in early- time spectra. A simila r 
result has already been reported bv lMatheson et all l)200ll l. 
whose data set is included here. However, also the variation 
of the FWHM from object to object increases in this di- 
rection, indicating particularly strong diversity in the ejecta 
geometry of BL-SNe Ic. An observed trend towards smaller 
FWHM at later epochs (see Fig. |8} could be explained by 
changes in the excitation conditions as a consequence of de- 
creasing densities, such as a transition to more local positron 
deposition as the dominant excitation mechanism. 



6 DISCUSSION OF INDIVIDUAL OBJECTS 

In the previous sections we have presented simple Gaussian 
fits to the [O i] AA6300, 6364 features in nebular spectra of 39 
stripped-envelope CC-SNe, and found substantial diversity 
in the line profiles with deviations from spherical symmetry 
in a majority of the objects. However, remarkable patterns 

2 Note that the BL-Ic SNe 1998bw and 2006aj, members of the 
narrow-core class, were discovered only after the detection of their 
associated GRBs. Given the particular geometry imposed by a 
GRB connection, a mild bias may be introduced into our statis- 
tics. Since we assert no claim to a strictly unbiased sample, we 
included these two objects in our analysis throughout this work. 
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Figure 8. FWHM of the one-component Gaussian fit to the [Ol] 
A6300 line as a function of time. Only selected SNe with good 
temporal coverage are displayed. 



can also be discerned within the zoo of line profiles. In the 
following the properties of some interesting individual SNe 
or groups of objects are discussed in more detail. 



6.1 SNe 19941, 1996N and 1996aq: 

Doppler-shifted blobs and neutron-star kicks? 

ISollerman et al.1 (|l99cf ) noticed a blueshift of [Ol] 
AA6300, 6364 in late-time spectra of SN 1996N. We confirm 
this result, obtaining an overall blueshift of 850 kms" 1 in 
one-component fits. The epochs of the spectra are too late 
180 d) to explain this with optical-depth effects. In a two- 
component setup, which yields a much better fit to the asym- 
metric line profile, the main component is nearly at rest, but 
the second one is blueshifted by ~3000 fans" 1 (Fig. [5] and 
Table El, prominent and broad (FWHM = 2300 kms" 1 , 
q~0.2 to 0.3). 

Here we stress the similarity of the [O i] line profiles 
of SNe 1996N and 19941 (although the latter was only ob- 
served at earlier phases, not later than 150 d past maxi- 
mum). Finding such an unusual [O i] profile in SN 19941 was 
unexpected, given that this is one of the most soundly stud- 
ied SNe Ic to date, and, to our knowledge, this fact has never 
been commented on in the literature. SN 1996aq features 
an apparently different, double-peaked [O i] line. However, 
the fitting suggests that the main difference is the width of 
the blueshifted component, which is significantly narrower 
in SN 1996aq (FWHM = 1000 kms" 1 ). This results in a sep- 
aration of the A6300 and A6364 lines and thus a two-horned 
appearance. 

The most convincing explanation for the observed line 
profiles in these three SNe is provided by blobs moving to- 
wards the observer at high velocity. If such a blob has an 
average composition and no enhanced excitation (as would 
be the result of an enhanced Co abundance), it has to carry 
substantial mass to account for the observed emission. In 



fact, in this simple scenario the mass fraction would be 
given by the fractional flux of the clump, a w . In SN 1996aq, 
a w is about 0.15, in SNe 19941 and 1996N 0.2-0.3. Mov- 
ing at a velocity > 3000 kms -1 , the blob carries enor- 
mous momentum, which has to be counter-balanced for 
the sake of momentum conservation. Since the remainder 
of the oxygen emission is centred at rest, the compensa- 
tion has to be provided by Fe- or Si-rich material, or the 
compact remnant of the core collapse. Strong neutron star 
kicks (up to several hundred kms -1 ) have inde ed been ob- 
served (e.g. ICordes. Romani fe Lundgrenl [l993h and repro- 
duced in simulations of anisotropic exp losions with dom- 
inant dipole (I = 1) mode in the ejecta llBurrows fe Havesl 
1 19961 ; IScheck et al]|2004 120061 ; iBurrows et al.ll2007f l. To esti- 
mate the kick velocities consistent with our measurements, 
we make a simple calcu lation for SN 1994 1. With a total 
ejecta mass of 1.2 M Q l|Sauer et all 120061 ) and an a w of 
0.22, the blueshifted blob contains 0.26 Mq if a homogeneous 
composition is assumed throughout the ejecta. To compen- 
sate the momentum, a typical neutron star of 1.3 Mq needs 
a kick velocity of ~600 kms -1 if moving along the line of 
sight, which is not an unreasonable number. 

6.2 SNe 1998bw and 2002ap: narrow cores? 

SNe 1998bw and 2002ap share a similar [O i] AA6300, 6364 
line profile with a narrow component on top of a much 
broader base. In SN 1998bw, this was attributed to emis- 
sion from a disk- or t orus-shaped oxygen distribution viewed 
nearl y from the top l)Mazzali et al.l[200ll ; iMaeda et al 1 l200l 
2006). Together with Fe lines being broader than [0 1] 
AA6300, 6364, this gave rise to the idea of a strongly aspheri- 
cal, jet-like e xplos ion. For the nebular spectra of SN 2002ap, 
Foley et al.l l)2005t ) remarked on a si milarity of the li ne pro - 
files with those of SN 1998bw, and IMazzali et ail (|2007aT ) 
suggested asphericity here as well, although the narrow peak 
might have been caused by a dense core in the ejecta. 

However, in both SNe the narrow components are 
redshifted with res pect to both the broad base and the 
rest wavelength Ao (jPatat et al.l l200ll ; iLeonard et al 1 120021 : 
iFolev et al.ll2003l ). From the best-fit parameters reported in 
Table ED mean redshifts of 586 ± 162 kms" 1 and 657 ± 62 
kms -1 with respect to Ao are inferred for SNe 1998bw and 
2002ap, respectively (cf. Fig|5J). Relative to the broad bases, 
mean offsets of 876 ± 129 kms" 1 and 500 ± 138 kms" 1 
are observed. Such large offsets (comparable to the narrow 
component's FWHM) are not expected for models with en- 
hanced central density or the bipolar-jet scenario favoured 
for SN 1998bw. Instead, the nebular [0 1] lines should be 
symmetric and centred at their rest wavelength unless the 
putative jet was rather one-sided. In that case, the lack of 
symmetry between the two hemispheres would explain the 
observed [0 1] line shifts. In this work, we conservatively 
classify SNe 1998bw and 2002ap as possible narrow cores 
(Fig. [2}, considering the presence of a blob another possibil- 
ity. Other SNe of the narrow-core group exhibit similar line 
offsets (cf. Table IA1|I , but mostly less pronounced than in 
SNe 1998bw and 2002ap. 

Inspired by the modelling of iTomita et ail l|2006h and 
IMazzali et all l|2007al ). we examined an alternative, possibly 
more physical configuration for SN 2002ap. This consisted 
of a jet (or, more generically, a bipolar explosion) viewed 
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Figure 9. Shifts of the broad line bases and the narrow cores of 
[O i] A6300 in SNc 1998bw and 2002ap as a function of time (cf. 
Table [XT> . 

strongly off-axis, but with an additional density enhance- 
ment at low velocity. To test the consistency with the ob- 
served [O i] line profile, we employed a three-component fit 
consisting of a strictly symmetric double peak plus a central 
narrow component (7 effective parameters; see Table |A2]) . 
and compared with the 'broad peak + narrow core' setup 
adopted for SN 2002ap throughout the rest of this paper. 

It turns out that the two configurations perform simi- 
larly well, the only advantage of the two-component fit being 
that one fewer parameter is involved. It is therefore difficult 
to decide on the basis of these fits which ejecta geometry (on- 
axis jet, spherical ejecta + dense core, spherical ejecta + blob 
or off-axis jet + dense core) is most likely. In the 'DP + NC' 
configuration, about 15 % of the emitting mass would be 
contained in the dense core. This is in good agreement with 
the size of the core that was added to reproduce the late-time 
light cu rve of SN 2002ap in the models (0.5 out of a total of 
3.0MM: lTomita et alJliooj ). The good quality of the fit and 
the consistency with sophisticated modelling make the 'DP 
+ NC scenario an attractive possibility for SN 2OO2ap0 

6.3 SNe 2003jd, 2000ew, 2004gt and 2006T: 
genuine double peaks and impostors 

Besides SNe 2003jd and 2006T, the prototypes of the double- 
peaked class, up to 5 other SNe have [O i] lines that agree 
with a DP configuration. The problem with most of these 

3 Also other SNe of the narrow-core class can be adequately fit 
with three components in a configuration similar as in SN 2002ap. 
This could help to solve the 'problem of missing double peaks', 
see Section [5~2|l . However, for the GRB-related SNe 1998bw and 
2006aj a 'DP + NC configuration — though providing a good fit 
(cf. Table IA2j l - is not expected to be correct, since these SNe 
are supposed to be viewed along the jet axis. Lacking modelling 
predictions for most other SNe, we have not explored the 'DP + 
NC option any further. 



objects is that the separation of the two fit components is 
smaller than in SNe 2003jd and 2006T, and often similar 
to the separation of the two [Ol] lines in the doublet. This 
leaves room for alternative interpretations, in particular the 
possibility that the two horns observed e.g. in SNe 2000ew 
and 2004gt may actually originate from the deblended A6300 
and A6364 lines of a single narrow, blueshifted component. 
In fact, in SN 2000ew the line profile is better reproduced by 
the latter configuration. This is the reason for its primary 
association with the AS class (see Fig. [5}. 

In SNe 2000ew and 2004gt also the intensity ratio of 
the two peaks appears inverted with respect to those of 
SNe 2003jd and 2006T, the blue peak being stronger than 
the red one. In a classical DP configuration (cf. Section [STTJ, 
however, the red peak will always be stronger since, for the 
given separation of the peaks, the A6364-line of the blue 
component blends with the A6300Tine of the red one. To 
circumvent this problem, a more complex ejecta structure 
with additional blueshifted emission on top of an otherwise 
symmetric profile may be assumed. This yields a good fit for 
SN 2004gt, cf. Table El Alternatively, the original toroidal 
oxygen distribution may be unchanged, but the redshifted 
emission component is damped owing to optically thick in- 
ner ejecta (cf. Section 0J. Since the spectra of SNe 2000ew 
and 2004gt are both relatively young (112 and 160 d, re- 
spectively) compared to those of SNe 2003jd and 2006T (cf. 
Table [TJ, this may indeed be a possibility. 

Due to its ejecta velocities and its double-peaked [O i] 
profile, SN 2003jd has been pro posed to be associat ed with a 
GRB viewed strongly off-axis (|Mazzali et al.ll2005l ). In that 
case, the 7- and X-ray emission of the GRB would not be 
seen because of the strong collimation of the jet. However, 
depending on the jet propagation model a radio afterglow 
would possibly be observable at late phases. SN 2003jd wa s 
not detected at radio wavelengths |Soderberg et al.1 12006). 
so that its association with a GRB remains uncertain. The 
second object with a similar [Oil line pro file, SN 2006T, was 
classified as SN lib l|Blondin et al.ll2006f ). This makes it per 
se a poor candidate for a GRB-SN, since the relativistic jet 
of a GRB would have to penetrate the He and H shells and 
probabl y die before reaching the surface. This suppor ts the 
view of iModiaz et all [|200Sh and iMaeda et all 0008) that 
strong asphericity is ubiquitous in core-collapse SNe, and 
not necessarily a signature of an association with a GRB 
(see also Section [5. 2p . 



7 THE PROFILE OF Mg I] A4571 

Hydrodynamic explosion models ( Maeda et al.l I2006T I sug- 
gest that Mg and O should have similar spatial distribu- 
tions within the SN ejecta, which may deviate significantly 
from those of heavi er elements such as Fe or Ca (see also 
iMazzali et al ] |2005h . This should result in the profiles of 
isolated Mg and O emission lines being simi lar, which so far 
has been shown t o hold for individual SNe (|Spyromiliolll994l ; 
iFolev et "aD I2003T 1 . Here we test this for our entire sample, 
examining the profile of the semi-forbidden Mg 1] A4571 line 
(3s 2 1 So - 3s3p 3 P°i) in all spectra with sufficient S/N in 
the wavelength range of interest. 

A direct comparison of the Mg 1] and [O 1] lines is hin- 
dered by the fact that, unlike Mgl], the [Ol] feature is a 
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Figure 10. Comparison of Mgl] A4571 and [Ol] AA6300, 6364 
line profiles (with a second component added to Mg i] artificially 
to account for the doublet nature of [Ol], see discussion). A sub- 
sample of objects with a good overlap is shown. The dot-dashed 
blue line is the modified Mg i] , the solid red line [O i] . 



19908 ■ A 




19901 






91 d r Jui 


V. 


90d 


i% 




y\ji\ 

l'v| 'V 






(■< ^ 






\ 'i*.J - 

V I;/ 


" 4- 






T ' 


\ I' 

A '' 


















.i / 











-0.02 0.00 0.02 -0.02 0.00 0.02 



v/c 



v/c 



1990W 
183d 


"ft , 


1997dq 
217d 


r» K. 
fi \\ 

\ \ ' 


2005N 
88d 




k 


i 


* w 

Juki's 

i 


1 


A \ 

\. \ 


/ 

ii ii 
A^ 






-0.02 


0.00 0.02 

v/c 


-0.02 


0.00 0.02 

v/c 


-0.02 


0.00 

v/c 


0.02 



2006gi 
148d 


AS* 


20071 A 
1 65d 




1 v. 
I \ i 


i V 


1 

. J'- 
- Wi;' 






-0.02 


0.00 0.02 


-0.02 0.00 0.02 



v/c 



v/c 



Figure 11. The same as Fig. 1101 but showing spectra with evident 
differences in the [Ol] AA6300, 6364 and the modified Mgl] A4571 
profiles. 



doublet. However, as described in Section [3.21 we assumed 
that in our nebular SN Ib/c spectra the A6300 and A6364 
lines have a ratio of 3 : 1. Therefore, to compensate we first 
isolated the Mgl] A4571 feature, subtracting a linearly fit 
background. Then we rescaled the Mgl] line to 1/3 of its 
initial intensity, shifted it by 46 A (equivalent to the 64 A 
offset of the two [O i] lines) and added it to the original pro- 
file. This modified Mgl] profile can then be compared with 
the observed [O i] feature. 

For most objects of our sample (~ 65%), even those with 
rather complex ejecta geometry, we find an impressive sim- 
ilarity of the Mg i] and [O i] line profiles within the noise 
level and the uncertainty in subtracting the background (see 
Fig. I10p . This indicates that the spatial distribution of Mg 
and O in the ejecta is generally similar. However, there are 
some noticeable exceptions with a poor match of the [O i] 
and modified Mg i] features. Examples are shown in Fig. 1111 
Most of these spectra are not very late, typically < 200 d, 
which leaves room for the following explanations: 

(i) In the affected objects, the O- and Mg-rich parts 
of the ejecta may indeed have different geometry owing to 
the chemical stratification of the progenitor star and the hy- 
drodynamics of the explosion. However, within this scenario 
differences in the line profiles should persist during the en- 
tire nebular phase. In SN 1998bw, for which a late nebular 
spectrum (376 d) is available, the differences visible at earlier 
epochs (< 200 d) are observed to vanish with time (Fig. ll2[) . 
SN 2007C undergoes a similar evolution even more rapidly. 

(ii) Alternatively, the Mg i] or the [O i] features may be 
contaminated by nebular emission lines of other elements. 
For the [Ol] feature the possibility of a contamination on 
the blue side has been discussed in Section [4] and found to 
be unlikely. Also, the fact that in most of the affected spectra 
the Mg i] line is broader than the [O i] line and changes more 
strongly with time (cf. Figs. [TT] and I12|l suggests that the 
Mg i] rather than the [O i] line may be blended with other 
lines at early epochs. 

(iii) The Mgl] A4571 line is located in a region shaped 
by many strong Fe II features during the photospheric phase. 
Hence, it is possible that the emission peak tentatively iden- 
tified as Mg i] is mostly produced by underlying photospheric 
Fe lines in some earlier spectra. This would not only explain 
the observed evolution in SN 1998bw (whe re Fe features 
are p articularly strong and persistent, see e.g. lMazzali et all 
2005), but also the better agreement of [Ol] and Mgl] in 
SN 2002ap early on: SN 2 002ap features the strongest Mgl] 
A4571 line ever observed (jFolev et al.ll2003h but only weak 
Fe, and it is hence not unexpected that Mgl] quickly domi- 
nates over photospheric residuals (Fig. I12|l . 



8 CONCLUSIONS 

We have studied the profiles of nebular emission lines in 
stripped-envelope CC-SNe, with the aim of constraining the 
ejecta morphology, and in particular the degree of aspheric- 
ity of the explosions. The study was based on 98 nebular 
spectra of 39 different SNe, some of which were not published 
before. The size of the sample gives it statistical significance. 
We have concentrated on the profile of [O i] AA6300, 6364, 
since this is usually the strongest feature in nebular SN Ib/c 
spectra, and not severely contaminated by other lines. We 
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Figure 12. Evolution of the [O i] AA6300, 6364 (solid red) 
and the modified Mgl] A4571 (dot-dashed blue) line profiles of 
SNe 1998bw, 2002ap and 2007C with phase. While in SN 2002ap 
the profiles are always similar, significant differences are apparent 
in the earlier spectra of SNe f 998bw and 2007C. 



decreasing phase. Geometrical effects and dust formation 
within the ejecta seem to be excluded as origins of the 
blueshift. Contamination from other elements may play a 
role in some, but not all SNe. Hence, residual opacity in the 
inner ejecta remains the most likely explanation for the ob- 
served shift, as photons emitted on the rear side of the SN 
are scattered or absorbed on their way through the ejecta, 
giving rise to a flux deficit in the redshifted part of spectral 
emission lines. The required opacity might be generated by 
a multitude of weak Fe transitions. 

A surprisingly good agreement of the profiles of [O i] 
AA6300, 6364 and Mgi] A4571 (modified to account for the 
doublet nature of the oxygen feature) was found in most SNe 
regardless of their morphological class. This indicates that 
the line profiles are indeed determined by the ejecta geome- 
try, and that Mg and O are similarly distributed within the 
SN ejecta. Deviations are mainly found in relatively early 
spectra, and we propose blending of the emerging nebular 
Mg i] emission with residual photospheric Fe n lines as a pos- 
sible reason for the differences. 



performed a multi-parameter Gauss-fitting of this feature in 
all our spectra, with the position, FWHM and intensity of 
the A6300 Gaussian being free parameters (and the A6364 
line duly added with fixed offset and intensity ratio of 1/3). 
Using this approach, a variable number of emission compo- 
nents, their widths and Doppler shifts could be identified. 
Compared to spectral modelling this method has the advan- 
tage of being fast, capable of dealing with complex profiles, 
and independent of an accurate flux calibration. 

Despite the large variety in line profiles encountered, 
we can divide the SNe into four morphologically different 
groups on the basis of the best-fit parameters: SNe with 
simple Gaussian line profiles (~ 18% of all objects), SNe 
with narrow line cores atop broad bases (~28%), objects 
with symmetric double peaks (~ 15%), and objects with ev- 
idence for blobs or overall asymmetric line profiles (~39%). 
Since this classification refers to the structure inferred from 
the fitting procedure, we believe it to be more relevant for 
the true ejecta geometry than a pure visual inspection of 
the line. The results of our analysis suggest that probably 
at least half of all SNe Ib/c are aspherical. The fraction of 
double peaks is too small for a jet to be a ubiquitous feature 
in a majority of the objects. Even among narrow-core SNe, 
only a small fraction of the objects may have a jet-like ejecta 
morphology. Instead, a central density enhancement appears 
to be a likely solution for many members of this class. A ma- 
jority of broad-line SNe Ic have narrow line cores, whereas 
most SNe lb exhibit asymmetric or multi-peaked line pro- 
files. 

Bulk shifts of the [O i] feature and strongly Doppler- 
shifted, massive blobs are observed in some objects, and ex- 
pected to be the signature of very one-sided explosions. If 
momentum conservation is provided by neutron-star kicks, 
the inferred kick velocities are compatible with those of the 
fastest-moving neutron-stars in the Galaxy. 

In spectra taken earlier than ~ 200 d after maximum 
light, a systematic blueshift of the [O i] AA6300, 6364 line 
centroids can be discerned, becoming more pronounced with 
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APPENDIX A: FIT PARAMETERS 

In the following tables the parameters of the best fits to 
the [O i] features of our spectra are reported. All numbers 
refer to the A6300 line, since the A6364 line is automati- 
cally added by the fitting algorithm at 1/3 of the stength. 
One- and two-component fits have been obtained for all 
SNe. They are listed in Table lA"T1 along with error estimates 
for the one-component fit parameters. Three-component fits 
have only been performed for a few objets. In SNe 1990B, 
20061d and 2007C the line profiles were too complex to be 
satisfactorily reproduced with fewer components, while for 
SNe 1997dq, 1998bw, 2002ap and 2006aj particular ejecta 
geometries were tested. The three-component fit parameters 
are reported in Table lA"2l 
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Table Al. Parameters™ of the best one-component (cols. 4—5) and two-component (cols. 6-12) fits of [Ol] A6300 - l at part. 



SN Date Sample region A FWHM \ 1 FWHM I ai A 2 FWHM 2 a 2 rel. RMS 6 
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6320 


5 


71. 


5 





55 


6270 


2 


72 


7 





45 





.840 


1997ef 


1998/09/21 


6166. 


1- 


-6470. 


1 


6294 


.9 


± 


5 


3 


120.2 ± 15.0 


6300 


6 


179. 


2 





72 


6290 


.8 


44 


.8 





.28 





.486 


1998bw 


1998/09/12 


6090. 


0- 


-6549. 


7 


6298 


.4 


± 


5 


.8 


185.5 ± 6.0 


6297. 


3 


191. 


2 


0. 


97 


6315 


1 


31 


.6 





.03 





.730 




1998/11/26 


6135. 


6- 


-6504. 


6 


6298 


.4 


± 


3 


.2 


146.1 ± 4.0 


6292 


8 


166. 


1 


0. 


87 


6315 


1 


42 


.6 





.13 





.124 




1999/04/12 


6134. 


6- 


-6500. 


3 


6299 


.3 


± 


3 


4 


120.6 ± 5.0 


6294 





142. 


8 


0. 


84 


6310 


.8 


39 


1 





.16 





.107 




1999/05/21 


6162. 


3 


-6460. 


7 


6297 


.6 


± 


3 


2 


109.8 ± 3.0 


6291 


7 


132. 


7 


0. 


82 


6308 


.2 


39 


.0 





.18 





103 


1999cn 


2000/08/04 


6198. 


2 


-6421. 


9 


6294. 


6 


± 


10 


.3 


108.4 ± 20.0 


6292 





123. 





0. 


96 


6303 


.7 


32 


.9 





04 





.785 


1999dn 


2000/09/01 


6214. 


1- 


-6455. 


2 


6309. 


6 


± 10 


.2 


94.3 ± 20.0 


6304 


.1 


112. 


9 


0. 


81 


6322 


.9 


33 


.8 





.19 





.920 
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Table Al. cont. Parameters" of the best one-component (cols. 4-5) and two-component (cols. 6-12) fits of [Ol] A6300 - 2 nd part. 



SN Date Sample region A FWHM Ai FWHMi oji A 2 FWHM 2 a 2 rel. RMS 1 

(1) (2) (3) (4) (5) (6) (7) (8) (9) (10) (11) (12) 



2000ew 


2001/03/17 


6189 


.4- 


-6420 


.2 


6284. 


3 


± 


3. 





99.3 ± 5.0 


6294. 


3 


109 


6 


0. 


79 


6266 


.3 


33 





0. 


21 





.496 


2002ap 


2002/08/06 


6082 


.7- 


-6529 


.4 


6297 


9 


± 


2. 


.3 


172.5 ± 5.0 


6296. 


(i 


179 


2 


0. 


96 


6312 


.3 


21 


9 


0. 


04 





.615 




2002/06/17 


6064 


.6- 


-6516 


.5 


6300 





± 


1 


.5 


174.0 ± 3.0 


6298. 


7 


181 


5 





95 


6312 


.5 


22 


5 


0. 


05 





.254 




2002/07/11 


6084 


.1- 


-6514 


.9 


6302 


4 


± 


1 


.5 


163.1 ± 1.0 


6301. 


6 


164 


6 


0. 


97 


6313 





20. 





0. 


03 





.324 




2002/08/09 


6103 


.9- 


-6503 


4 


6304 


6 


± 


1 


.5 


157.5 ± 1.0 


6303. 


5 


165 


7 


0. 


96 


6314 


.1 


21 


6 


0. 


04 





.164 




2002/10/01 


6094 


.6- 


-6513 


3 


6305 


.4 


± 


1 


.5 


150.4 ± 2.0 


6304. 


1 


157 


9 


0. 


95 


6314 


.4 


2(1 


2 


0. 


05 





.216 




2002/10/09 


6104 


.7- 


-6511 


.7 


6305 





± 


1 


.5 


152.2 ± 4.0 


6303. 


9 


160 


6 





95 


6313 


.5 


20. 


3 


0. 


05 





.366 




2002/10/14 


6104 


.0- 


-6511 


.2 


6306 


.4 


± 


1 


.5 


153.9 ± 3.0 


6304. 


2 


166 


4 





93 


6316 


.2 


31. 





0. 


07 





.435 




2002/11/06 


6149 


.5- 


-6500 


.5 


6305 


.7 


± 


1 


.5 


149.0 ± 1.0 


6304. 


7 


157 


8 


0. 


95 


6312 


.5 


22 


2 


0. 


05 





.247 




2003/01/07 


6137 


.7- 


-6519 





6306 


9 


± 


2. 


.3 


143.0 ± 2.0 


6305. 


9 


150 


8 





95 


6313 


.7 


21 


7 


0. 


05 





.373 




2003/02/27 


6145 


0- 


-6490 





6309 


.5 


± 


1 


.5 


140.0 ± 2.0 


6309. 





140 


8 


0. 


97 


6315 


.4 


17. 


1 


0. 


03 





.534 


2003jd 


2004/09/11 


6087 


.3- 


-6505 


.8 


6299 


1 


± 


3. 


.6 


219.6 ± 8.0 


6232. 





116 


7 


0. 


50 


6357 


2 


87. 





0. 


50 





.435 




2004/10/18 


6150 


.7- 


-6488 


.2 


6304 


.4 


± 


6. 


3 


187.5 ± 15.0 


6251. 


4 


87. 








52 


6358 


.5 


67. 


8 


0. 


48 





.584 


2004aw 


2004/11/14 


6124 


.1- 


-6498 


.6 


6296 


.1 


± 


3. 


.3 


118.7 ± 7.0 


6293. 


1 


133 





0. 


89 


6307 


.6 


37. 





0. 


11 





.507 




2004/12/08 


6191 


.6- 


-6430 


.1 


6289 


.2 


± 


(i 


.6 


91.0 ± 14.0 


6287. 


5 


94. 


7 





94 


6302 


.9 


34 


8 


0. 


06 





.882 




2005/05/11 


6199 


.0- 


-6442 





6298 





± 


6. 


.6 


91.4 ± 16.0 


6294. 





101 


8 





87 


6310 


.8 


32. 


7 


0. 


13 





.651 


2004gt 


2005/05/24 


6174 


.8- 


-6485 


.3 


6285 


.2 


± 


3. 


1 


125.1 ± 3.0 


6260. 


7 


72. 


5 


0. 


60 


6335 


.8 


87. 


1 


0. 


40 





.490 


2005N 


2005/01/21 


6160 


0- 


-6476 


.3 


6268 


(i 


± 


6. 


.5 


117.3 ± 7.5 


6284. 





126 


3 





77 


6241 


.2 


46. 


8 


0. 


23 





.353 


2006F 


2006/11/16 


6223 


,0- 


-6450 


.9 


6296 


.9 


± 


6. 


3 


75.8 ± 20.0 


6314. 


6 


84 


7 


0. 


64 


6281 


,4 


38. 


9 


0. 


36 





.508 


2006T 


2007/02/18 


6203 


.0- 


-6452 


.4 


6300 


.5 


± 


8. 


.0 


123.6 ± 16.0 


6266. 


8 


50. 


7 





51 


6338 





57. 


9 


0. 


49 





.306 


2006aj 


2006/09/19 


6101 


0- 


-6498 


.7 


6294 


(i 


± 


6. 


.1 


188.5 ± 9.0 


6295. 





192 


2 





99 


6301 


.8 


IK. 


1 


0. 


01 





.891 




2006/11/27 


6096 


.0- 


-6515 


.7 


6303 


.5 


± 


(i 


.7 


203.4 ± 11.0 


6304. 


5 


222 


4 


0. 


93 


6300 


.8 


35. 


3 


0. 


07 





.588 




2006/12/19 


6119 


.2- 


-6507 





6305 


.2 


± 


7. 


.7 


177.1 ± 12.0 


6306. 


2 


197 


4 





92 


6302 


.3 


22 


4 


0. 


08 





.410 


2006gi 


2007/02/10 


6138 


.3- 


-6433 





6278. 


.5 


± 


4. 





107.0 ± 2.0 


6280. 


1 


114 


5 





93 


6267 


.7 


35. 


2 


0. 


07 





.458 


20061d 


2007/07/17 


6176 


.4- 


-6492 


.5 


6285 


.5 


± 


2 


.2 


97.2 ± 5.5 


6304. 


8 


104 





0. 


70 


6259 


.7 


43 


9 


0. 


30 





.572 




2007/08/06 


6173 


.5- 


-6491 


.2 


6287. 


.6 


± 


3. 


.1 


102.1 ± 9.0 


6305. 


7 


115 


1 


0. 


72 


6263 


.5 


43. 


8 


0. 


28 





531 




2007/08/20 


6175 


.5- 


-6482 


.4 


6290 


.3 


± 


2. 


6 


92.7 ± 7.0 


6296. 


3 


87. 


3 





89 


6251 


.9 


33 


9 


0. 


11 





936 


2007C 


2007/05/17 


6131 


0- 


-6501 


.4 


6278. 


.2 


± 


3. 


6 


122.5 ± 2.0 


6301. 


1 


92 


3 


0. 


69 


6235 





52 


4 


0. 


31 





.446 




2007/06/20 


6134 


0- 


-6502 


.5 


6285 


3 


± 


2 


.6 


121.5 ± 1.5 


6299. 


6 


103 


2 





SI 


6237 


.6 


41 


6 


0. 


19 





.344 


20071 


2007/06/18 


6157 


.1- 


-6494 


.7 


6300 


.3 


± 


5. 


.3 


118.8 ± 1.0 


6305. 


3 


126 


1 





90 


6277 


.8 


35. 


8 


0. 


10 





.122 




2007/07/15 


6159 


.2- 


-6504 


.8 


6303 


.4 


± 


5. 


.3 


1X1.9 ± 1.5 


6308. 


2 


116 


3 





90 


6278 


.8 


33. 


7 


0. 


10 





129 



a Central wavelengths (Aj), full widths at half maximum (FWHMj) and relative strengths (Qj) of the fit components (cf. Section l3.2| . The sample region, A 
and FWHM arc given in units of A. 

b RMS of the fit, relative to the RMS of the one-component fit. 

c The fitting was performed in Ho:-subtractcd spectra. For this purpose, a smoothed symmetric Ha profile was subtracted, constructed by reflecting the red 
wing at the rest wavelength (see also[Patat ct al. 1995T) ■ 



Table A2. Parameters' 2 of the best three-component fits of [Ol] A6300 for selected spectra with complex line profiles or particular 
geometric configurations. 



SN 


Date 


Sample region 


Ai 




FWHMi 


"1 


A 2 




FWHM 2 


c*2 


A 3 


FWHM 3 


a 3 


RMS 6 


Config. 


(1) 


(2) 


(3) 


(4) 




(5) 


(6) 


(7) 




(8) 


0) 


(10) 


(11) 


(12) 


(13) 


(14) 


1990B 


1990/04/18 


6132.1-6499.7 


6319 


.9 


180.7 


0.71 


6235 


.5 


50.0 


0.26 


6303.1 


13.0 


0.03 


0.140 


clumpy 


1997dq 


1998/05/30 


6205.0-6446.5 


6272 


.4 


63.3 


0.49 


6326 


.9 


63.3 


0.49 


6298.4 


12.3 


0.02 


0.480 


DP + NC 


1998bw 


1999/05/21 


6162.3-6460.7 


6259 


.4 


118.2 


0.38 


6319 


.4 


118.2 


0.38 


6307.8 


45.0 


0.24 


0.175 


DP + NC 


2002ap 


2002/10/01 


6094.6-6513.3 


6251 


.8 


102.6 


0.42 


6352 


.1 


102.6 


0.42 


6313.5 


42.3 


0.15 


0.170 


DP + NC 


2004gt 


2005/05/24 


6174.8-6485.3 


6243 


.8 


63.7 


0.37 


6338 


6 


63.7 


0.37 


6279.3 


44.3 


0.26 


0.229 


DP + blob 


2006aj 


2006/12/19 


6119.2-6507.0 


6256 


.9 


151.4 


0.45 


6353 


.3 


151.4 


0.45 


6302.2 


24.1 


0.09 


0.402 


DP + NC 


20061d 


2007/07/17 


6176.4-6492.5 


6320 


.8 


121.4 


0.48 


6261 


.6 


40.1 


0.40 


6301.6 


20.0 


0.12 


0.161 


clumpy 




2007/08/06 


6173.5-6491.2 


6320 


.4 


136.6 


0.53 


6264 


.2 


39.4 


0.36 


6302.1 


21.1 


0.11 


0.212 


clumpy 




2007/08/20 


6175.5-6482.4 


6322 


3 


117.5 


0.47 


6264 


.4 


46.5 


0.37 


6303.2 


19.9 


0.16 


0.343 


clumpy 


2007C 


2007/05/17 


6131.0-6501.4 


6299 





127.1 


0.56 


6239 


.5 


54.2 


0.29 


6301.1 


43.1 


0.15 


0.088 


NC + blob 




2007/06/20 


6134.0-6502.5 


6297 


.7 


125.6 


0.68 


6242 





43.5 


0.19 


6304.4 


41.4 


0.13 


0.060 


NC + blob 



a Central wavelengths (A^), full widths at half maximum (FWHMi) and relative strengths (a^) of the fit components (cf. Scction l3.2| . The sample region, A 
and FWHM arc given in units of A. 

b RMS of the fit, relative to the RMS of the one-component fit (cf. Table |~\T1 . 
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